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We propose a new class of gravity-matter theories, describing R -h R^ gravity interact¬ 
ing with a nonstandard nonlinear gauge field system and a scalar “dilaton”, formulated 
in terms of two different non-Riemannian volume-forms (generally covariant integration 
measure densities) on the underlying spacetime manifold, which are independent of the 
Riemannian metric. The nonlinear gauge field system contains a square-root \/ —F^ of 
the standard Maxwell Lagrangian which is known to describe charge confinement in 
fiat spacetime. The initial new gravity-matter model is invariant under global Weyl- 
scale symmetry which undergoes a spontaneous breakdown upon integration of the non- 
Riemannian volume-form degrees of freedom. In the physical Einstein frame we obtain 
an effective matter-gauge-field Lagrangian of “k-essence” type with quadratic depen¬ 
dence on the scalar “dilaton” field kinetic term X, with a remarkable effective scalar 
potential possessing two infinitely large fiat regions as well as with nontrivial effective 
gauge coupling constants running with the “dilaton” (p. Corresponding to the each of the 
two fiat regions we find “vacuum” configurations of the following types: (i) p = const 
and a non-zero gauge field vacuum \/ —F^ ^ 0, which corresponds to a charge confining 
phase; (ii) X = const (“kinetic vacuum”) and ordinary gauge field vacuum V—= 0 
which supports confinement-free charge dynamics. In one of the fiat regions of the ef¬ 
fective scalar potential we also find: (iii) X = const (“kinetic vacuum”) and a non-zero 
gauge field vacuum V—^ 0, which again corresponds to a charge confining phase. 
In all three cases the spacetime metric is de Sitter or Schwarzschild-de Sitter. Both “ki¬ 
netic vacuums” (ii) and (iii) can exist only within a finite-volume space region below a 
de Sitter horizon. Extension to the whole space requires matching the latter with the 
exterior region with a nonstandard Reissner-Nordstrom-de Sitter geometry carrying an 
additional constant radial background electric field. As a result we obtain two classes of 
gravitational bag-like configurations with properties, which on one hand partially paral¬ 
lel some of the properties of the solitonic “constituent quark” model and, on the other 
hand, partially mimic some of the properties of MIT bags in QCD phenomenology. 

Keywords: modified gravity theories, non-Riemannian volume-forms; global Weyl-scale 
symmetry spontaneous breakdown; fiat regions of scalar potential; charge confining non- 
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linear gauge field system; gravitational bags 

PACS numbers: 04.50.Kd, ll.30.Qc, 98.80.Bp, 95.36.+X 

1. Introduction 

To understand the basic features of hadronic physics we are unavoidably lead to 
the existence of different phases of the gauge theory. Furthermore, these different 
phases are apparently associated with different values of th e vac uum energy density, 
as described for example in the MIT bag model of hadrons 4121 Inside the MIT bag 
quarks and gluons propagate almost freely and there the vacuum energy density 
is big defining the so called bag constant.Outside the bag there is no propagation 
of either quarks or gluons, and in the MIT bag model the outside vacuum energy 
density is set to zero (any choice for the outside energy density is possible, if we 
ignore gravity, since while ignoring gravity only the difference of the energy den¬ 
sities inside and outside the bag is of significance). For a vacuum state p = —p, 
where the vacuum pressure inside the bag is negative while zero outside, an empty 
bag therefore tends to implode. When the bag is filled with particles, the positive 
pressure of the particles stabilizes the bubble at a certain radius. 

It is interesting that a similar “two phase” structure defined via two vastly 
different scales of the vacuum energy density does also appear in cosmology. Indeed, 
according to the most accepted scenario of the early universe - the inflation picture 
there was at the beginning a large vacuum energy density. Together wi th t his, 
for a description of the present slowly accelerated phase of the universe EHIS] Qjjg 
employs a small vacuum energy density. 

A scenario of continuously connecting an inflationary phase to a slowly accelerat¬ 
ing universe through the evolution of a single scalar field - the so called “quintessen¬ 
tial inflation” scenario - has been first proposed in Ref. 1141 whic h tr iggered active 
further development (models based on generalized F{R) gravity based on the 
k-essence framework - see Ref. [20l based on the “variable gravity” mod 

and containing an extensive list of references to earlier work on the topic - see 
Ref. [22l see also Ref. 1^ . 

In the cosmological context we have been able to construct models providing 
a unified scenario where both an inflation and a slowly accelerated phase for the 
universe can appear naturally from the existence of two infinitely large flat regions 
in the effective scalar field potential with vastly differ ent scales which we derive sys¬ 
tematically from a Lagrangian action principle l ^4 | 25 | have constructed 

a new kind of globally Weyl-scale invariant gravity-matter action within the first- 
order (Palatini) approach formulated in terms of two different non-Riemannian 
volume-forms (generally covariant integration measure densities on the pertinent 
spacetime manifold independent of the Riemannian metric). The principal feature 
is the requirement of global Weyl-scale invariance and the choice of different scaling 
properties of the two non-Riemannian measures (volume elements) which dictates 
the precise form of the terms in the action. In this new theory there is a single 
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scalar field with kinetic terms coupled to both non-Riemannian measures, whereas 
the standard Einstein-Hilbert term R and the i?^-term couple each to a different 
non-Riemannian measure. 

Global Weyl-scale invariance is spontaneously broken upon solving part of the 
equations of motion corresponding to the auxiliary antisymmetric tensor gauge fields 
of maximal rank defining the two non-Riemannian measures - due to the appearance 
of two arbitrary dimensionful integration constants. The l atter p roduce a remark¬ 
able effect on the resul ting physical Einstein-frame theory - we find there an 

effective k-essence LLSHiMJ type of theory, where the effective scalar field potential 
has two infinitely large flat regions corresponding to the two accelerating phases 
of the universe - one for large negative values of the scalar field with a very large 
height corresponding to the early universe, and another one for large positive values 
of the scalar field with a very low height corresponding to the universe of the present 
epoch. 

Since the construction was based on geometrical and symmetry considerations, 
which are very general, one may think that a similar model can be constructed for a 
different physical application, f.e, the phases of a gauge theory. As we will see here, 
it is indeed possible to obtain a phase structure of confinement and deconfinement 
related to what the MIT bag model suggests. 

To this end we will couple the above new type of gravity-matter theory defined 
in terms of the two different non-Riemannian measures and containing Rf gravity 
term to a special kind of nonstandard nonlinear gauge field model (for the analogous 
situation in the less general case of gravity-matter models with one non-Riemannian 
and one standard Riemannian integration measures, see Ref. 1261) . 

Namely, let us consider Abelian gauge fields whose Lagrangian contains both 
the standard Maxwell Lagrangian as well as the nonstandard 

square-root of the latter \/—F^) coupled to the two different non-Riemannian 
measures in a globally Weyl-scale invariant form. In flat s pac e-time the \/—F^- 
term is known to describe dynamics of charge confinement 12Z1-131J the 

non-zero vacuum va lue of V—F^. The latter is an explicit realization of an earlier 
proposal by ‘t Hooft who argued that the energy density of electrostatic field 
configurations in the low-energy description of confining quantum gauge theories 
must be a linear function of the electric displacement field in the infrared region 
(the latter appearing as a quantum “infrared counterterm”). Below, in Appendix B 
we extend the flat spacetime proof in Ref. [351 about the charge confining property 
of the V—F^-term to the case of curved static spherically-symmetric spacetimes. 

For further interesting properties of gravity-matter theories involving the 
“square-root” Maxwell term V—F^ (black holes with confining electric potential, 
new mechanism for dynamical generation of cosmological constant, charge-hiding 
and charge-confining via “tube-like” wormholes), see Refs. [34ll36l 

Let us remark that one co uld start with th e non-Abelian version of the nonlin¬ 
ear gauge field system with ■ Since we will be interested in static 

spherically symmetric solutions, the non-Abelian theory effectively reduces to an 
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Abelian one as pointed out in Refs. [27ll^ 

In the present context after including the coupling of the non-Riemannian- 
measures-modified gravity-matter theory to the nonlinear gauge field with \/— 
in the pertinent physical Einstein frame we obtain an effective matter Lagrangian 
again of “k-essence” type with quadratic dependence on the ip kinetic term X of 
the scalar “dilaton” field with a remarkable effective scalar potential possessing 
two infinite flat regions with different energy scales. In addition we get a nontriv¬ 
ial coupling of the nonlinear gauge field to the “dilaton” kinetic term 
All terms are multiplied by nontrivial “dilaton”-dependent coefficient functions, in¬ 
cluding nontrivial effective gauge coupling constants running with Lp. An important 
observation is their “flatness” (constancy w.r.t. running (p) in both infinitely large 
flat regions of the effective potential. 

We study the static spherically symmetric “vacuum” configurations correspond¬ 
ing to each of the two flat regions. In all cases the gravitational part is de Sitter type 
(de Sitter or Schwarzschild-de Sitter) with effective cosmological constant whose 
value is determined by the height of the total effective “dilaton” potential in each of 
the flat regions. The latter includes, apart from the purely scalar “dilatonic” ones, 
also additional contributions due to (possible) non-zero vacuum values of \/— 
and X. 

The static spherically symmetric “vacuum” configurations of the “dilaton” ip 
and the nonlinear gauge field are of the following types: 

(i) p = const, i.e., X = 0 and a non-zero gauge field vacuum V—F^ ^ 0, the 
latter corresponding to a confining phase - these solutions exist in both flat regions 
of the effective scalar potential (the one for large negative values of p and the other 
one for large positive values of p) 

(ii) X = const (“kinetic vacuum”) and ordinary gauge field vacuum V—F^ = 0, 
which supports confinement-free charge dynamics - this solution exists both in the 
flat region of the effective scalar potential for large positive values of p as well as 
for a special value of one of the scalar potential’s parameters also in the flat region 
for large negative values of p. 

(iii) X = const (“kinetic vacuum”) and a non-zero gauge field vacuum \/—F"^ ^ 
0, which again corresponds to a confining phase - this solution exist in the flat 
region of the effective scalar potential for large negative values of p for generic 
scalar potential’s parameter values. 

An important point here is that both “kinetic vacuums” (ii) and (iii) do not 
represent themselves genuine vacuum configurations, since they are defined only 
within a finite-volume space region below the de Sitter horizon. In order to obtain 
a well-defined static spherically symmetric configuration over the whole spacetime 
we need to match at the de Sitter horizon the “kinetic vacuums” (ii) and (iii) 
living in the interior de Sitter region to the exterior region with a nonstandard 
Reissner-Nordstrom-de Sitter geometry which carries an additional constant radial 
background electric field. Studying the vacuum energy densities inside and outside 
the de Sitter horizon shows that the inside energy density is higher than the outside 
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one. Thus, the fully extended to the whole spacetime “kinetic vacuums” (ii) and 
(iii) represent gravitational bag-like configurations where: 

(a) Th e type (ii) gravitational bag mimics some of the properties of the MIT bag 
model finite volume space region with deconfinement and large energy density 
versus infinite volume exterior region with confinement and low energy density; 

(b) Both type (ii) and type (iii) gravitational bag s resemble some of the prop¬ 
erties of the solitonic “constituent quark” model 123 - they are charged and carry 
“color” flux to infinity. 

The plan of the paper is as follows. In the next Section 2 we describe in some 
detail the general formalism for the new class of gravity-matter systems defined in 
terms of two independent non-Riemannian integration measures. In Sections 3 and 
4 we describe the properties of the two flat regions of the Einstein-frame effective 
scalar potential and derive the static spherically symmetric vacuum configurations. 
In Section 5 we construct static spherically symmetric solutions representing gravi¬ 
tational bag-like configurations. We conclude in Section 6 with some discussions. In 
Appendix A we briefly outline the canonical Hamiltonian treatment of the modified 
gravity-matter models with two non-Riemannian spacetime volume-forms, which 
elucidates the physical meaning of the auxiliary fields defining the non-Riemannian 
volume-forms. In Appendix B following Ref. [28] we show that the presence of the 
“square-root” Maxwell term \/—generates confining effective potential between 
quantized charged fermions in static spherically symmetric spacetimes. 


2. Gravity-Matter System Coupled to Charge-Confining Nonlinear 
Gauge Field — A Formalism With Two Independent 
Non-Riemannian Volume-Forms 


We shall consider the following nonstandard gravity/nonlinear-gauge-field/matter 
system with an action of the general form involving two independent non- 
Riemannian integration measure densities generalizing the models studied in 
Refs. 113 mi (for simplicity we will use units where the Newton constant is taken as 
CNewton = I/IGtt): 


S= / d‘^x<^>i{A) 




+ / d^X^2iB) 


+ £i?2 


m) 




( 1 ) 


Here the following notations are used: 


• $i(A) and $2(5) are two independent non-Riemannian volume-forms, i.e., 
generally covariant integration measure densities on the underlying space- 
time manifold: 

$i(A) = , $2(5) = , (2) 

defined in terms of field-strengths of two auxiliary 3-index antisymmetric 
tensor gauge fields. $ 1,2 take over the role of the standard Riemannian 
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integration measure density ^/—g = ^/--detYg^^ in terms of the space- 
time metric 

• R = g^''Rfi^{T) and Rf^^{T) are the scalar curvature and the Ricci tensor 
in the first-order (Palatini) formalism, where the affine connection is a 
priori independent of the metric g^i,. Note that in the second action term 
we have added a R^ gravity term (again in the Palatini form). Let us recall 
that R + R? gravity within the second order formalism (which was also the 
first inflationary model) was originally proposed in Ref. |38l 

• denote two different Lagrangians of a single scalar matter field (“dila- 
ton”) and of a Abelian gauge field potential of the form: 

- y\/-F2(g) , V(ip) = fiexp{-aip}{?) 

= _!ie-<^Vgi^''d^ipd^ip + U{ip) - ^F^{g) , U{(p) = f2 e^p{-2a(p}{^) 

where: 

F^ig)=Ff,,F,xgf^^g''^ , F^, = . (5) 

Here q;,/i ,/2 are dimensionful positive parameter, whereas 6 is a dimen¬ 
sionless one. The choice of the scalar potentials in ©-(HI) is similar to the 
choice in Ref. [39l 

• d>(77) indicate the dual field strength of a third auxiliary 3-index antisym¬ 
metric tensor gauge field: 

HH) = , ( 6 ) 

whose presence is crucial for non-triviality of the model. 

Concerning the explicit form of the non-Riemannian integration measure densi¬ 
ties (El let us note that any of the pertinent auxiliary 3-index antisymmetric tensor 
gauge helds, for instance, A^i,x can be in particular parametrized in terms of 4 
auxiliary scalar fields 

Af,^\ = ^eijKL (j)^dx(j)^ , (7) 

so that: 

^’i(A) = dxij}^ = det ll|^ll (8) 

acquires the form of a Jacobian. In a recent study of general relativity as an 
extended canonical gauge theory a similar Jacobian representation of the covariant 
integration measure has appeared in terms of additional scalar fields. However, 
unlike the present case in the construction of Ref. [40] the additional scalar fields 
enter also in the proper Lagrangian. 

In what follows we will stick to the representation EJ. 
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The scalar field potentials and the separate locations of the standard Maxwell 
and the square-root Maxwell gauge field terms have been chosen in such a way that 
the original action o is invariant under global Weyl-scale transformations: 


t ’ T’ + ~ In A , —>• , 






(9) 


The equations of motion resulting from the action o are as follows. Variation 
of H]) w.r.t. affine connection 

Jd^x + 2e ^ i?) - v^dr;:,) = o (lo) 

gives, following the analogous derivation in the Ref. [39l that becomes a Levi- 
Chevitta connection: 

K\ = Kxia) = ^g^'" {dug\K + dxg^K - d^g^x) , (H) 

w.r.t. to the Weyl-rescaled metric 


g^I.v = (xi + 26X277)5^1^ 


, Xi = 


$i(A) 

V^g 


X2 


W) 
V^g ■ 


( 12 ) 


Variation of the action © w.r.t. auxiliary tensor gauge fields A^^x, 7?^i/A and 
yields the equations: 




= 0 


+ ei?2 


$(iJ) 




= 0 


whose solutions read: 



(13) 


= y 2 = const , R+lA'^ = —Ml = const , L^^^-|-eJ7^+= —M 2 = const . 

(14) 

Here Mi and M 2 are arbitrary dimensionful and X 2 arbitrary dimensionless in¬ 
tegration constants. The appearance of Mi, M 2 signifies dynamical spontaneous 
breakdown of global Weyl-scale invariance under ([9|) due to the scale non-invariant 
solutions (second and third ones) in (|T4|l . 

It is also very instructive to elucidate the physical meaning of the three arbitrary 
integration constants Mi, M 2 , X 2 from the point of view of the canonical Hamilto¬ 
nian formalism. Namely, as shown in the Appendix A Mi, M 2 , X 2 are identified as 
conserved Dirac-constrained canonical momenta conjugated to (certain components 
of) the auxiliary maximal rank antisymmetric tensor gauge fields A^^x, B^,jx,H^vX 
entering the original non-Riemannian volume-form action ©. 

Varying o w.r.t. g^i, and using relations (|14l) we have: 


Xi 


Rn 


+ i - T«)] - 1x2 [Tl^+g>.. {eR^ + M 2 ) - 2RR^,] = 0 , (15) 
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f -j 2^ 

where xi and X 2 are defined in (Ha, and T/iiy’ ' are the energy-momentum tensors 
of the scalar-|-gauge field Lagrangians with the standard definitions: 




(16) 


Taking the trace of Eqs. m and using again second relation (imi we solve for 
the scale factor xi- 

T(2)/4 -t M 2 


Xi _ Ti^)/2 - Ml " 

where 

Using second relation m Eas. (fT51) can be put in the Einstein-like form: 

' \tI^J+ g^jM2 + e{L^^^+Mif 


(17) 


2 xiUL 


where: 


rj = 1 - 2e -k Ml) . 


(18) 


(19) 


Let us note that m, upon taking into account second relation m and m , can 
be written as: 

9ilv =Xi^9^iu ■ ( 20 ) 

Now, we can bring Eas. lfT^ into the standard form of Einstein equations for the 
rescaled metric g^i, ([20)l . i.e., the Einstein-frame equations: 

R^^M - \9^..R{9) = \Tf. ( 21 ) 

with energy-momentum tensor corresponding according to the definition (1161) : 


Tf, = g^,L,s - 2 


d 


dgt^'' 


-^eff 


to the following effective (Einstein-frame) scalar field Lagrangian: 

Ueff = ^ } ■ 


( 22 ) 


(23) 


XiU L 

In order to explicitly write Leff in terms of the Einstein-frame metric 5 ^ 1 / ()20l) 
we use the short-hand notation for the scalar kinetic term: 




(24) 


and represent in the form: 

=xiUX-U-xiU^y^F^ , =xinhe-<^^X+U-{xi^f A^F\g) , 

(25) 
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with V and U as in ©-O- 

From Eqs. (El) and (HH), taking into account (ESI), we find: 


(F - Ml) 




2 X 2 


[/ + M 2 + e(F-Mi )2 


1 - X 2 ( 


F- Ml 


-2e)X-ex2fo\/^FW) 


Upon substituting expression EE) into 
Einstein-frame matter Lagrangian: 


(26) 

we arrive at the explicit form for the 


Leff = Aiip)X+B{ip)X^-Uesiip)- 


4eeff(v5) 


\/-F'^{g)-eX 2 foA{tp)X^J-F'^ig) . 


The coefficient functions in (l?7l) read: 

Ai^) = l-AU,si^) - ® 


eX 2 - 


2(F((^)-Mi)J 
whereas the effective scalar field potential reads: 


Biip) = ex2-^Ues{<f) 


eX 2 - 


(27) 

X 26 e-“^ i 2 


2(F(^) - Ml) 
(28) 


(F-Mi )2 - Ml)' 

4X2 

U-k M 2 -f e(F - Mi)2 

CM 

1 

/ 2 e- 2 “v + M 2 + e(/ie-“'^ - Mi )2 

(29) 


where the explicit form of F and U are inserted. Eurther, the original gauge 

coupling constants are here replaced by (/j-dependent effective coupling constants: 

/eff(7>) = /o(l - 4eX2t/eff(V5)) 

/ 2 e- 2 “‘^ -p M 2 


= /o 

1 


/2e-2av + M 2 + e{fie-°“P - Mi)2 ’ 
1 


^eff 


iv) 


= X2 


+ e/o(l-4ex2C/eff(v?)) 


= X2 


-y + e/o 


f2e-^°^^ + M 2 


Le^ 


(30) 


(31) 


'°/2e-2“v5 + M2-f e(/ie-“‘^-Mi)2j ' 

Let us recall that the dimensionless integration constant X 2 systematically appearing 
in most relations is the ratio of the original second non-Riemannian integration 
measure to the standard Riemannian one El- 

We observe that even if we start with no standard Maxwell kinetic term for the 
gauge field, i.e., taking the limit —>■ oo in the original action 0 - 0 , we neverthe¬ 
less obtain a dynamically induced Maxwell term in the Einstein-frame action (12711 
with effective running charge according to EB: 

1 (/ 2 e- 2 “‘^ + M 2 ) 


xn-x = ^X2fo- 


(32) 


[/2e-2“v= + M 2 + e(/ie-“^ - Mi)2] ’ 

Erom (15^ we see that dynamical Maxwell term generation is a cumulative effect of 
the simultaneous presence of the “confining” gauge field term y/—F'^ and the 
gravity term. 
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Fig. 1. Qualitative shape of the effective scalar potential Ug{i{ip) i29t for Mi > 0. 


3. Flat Regions of the Effective Scalar Potential and Nontrivial 
“Vacuum” Solutions 

The explicit expressions for the effective potential Ueff{ip) (1^51) . the scalar (“k- 
essence”) kinetic terms’ coefficient functions A{lp) and R(v?) (1^51) and the effec¬ 
tive gauge coupling constants (ISH - d^Tl) reveal the following crucial feature of the 
Einstein-frame matter Lagrangian (Elll-dSIl): the presence of two infinitely large flat 
regions - one for large negative and one for large positive values of the scalar field 
ip, where all of the above are essentially constant w.r.t. ip. 

Depending on the sign of the integration constant Mi we obtain two types of 
shapes for the effective scalar potential Ueff{ip) (l29l) depicted on Fig.l. and Fig.2. 

For large negative values of ip we have for the effective potential and the coeffi¬ 
cient functions in the Einstein-frame matter Lagrangian (EZll-dSIl): 


Ues{ip) ^ [/(-) = 


/lV/2 


A{ip) ~ A(_) = 


1 + Wl/f2 


eisi^P) - e(-) = 


l + e/iV/2 

l + e/iV/2 


B{ip) ~ R(_) = ex2 


4X2(1 + e/i 7/2) 

1 + 6 / 1 7/2 - 6V4e/2 

l + e/iV/2 


X 2 1 + e/iV /2 + e^e/^ 


, fes{(p) ^ /(-) = 


/o 


l + e/iV /2 


, (33) 
, (34) 
• (35) 


This will be called “(-) flat region”. In the second flat region for large positive ip, 
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^eff 



Fig. 2. Qualitative shape of the effective scalar potential I129II for Mi < 0. 


which will be called “(+) flat region” we have: 

Ues{(p) ^ 17(+) = 


M2/M2 


A{^) ~ A(+) = 


M 9 , 


elsi^P) - e(+) = 


B{(p) ~ B(+) = ex2 
fes{>f) ^ /(+) = 


4x2(1 + eMl/M2) 
M 2 


M 2 + eMf 
1 + eMi jM-i 

Vi 1 + eM2/M2 + e2e/2 ’ W “ i + 


M 2 + eMf 

fo 


,(36) 

,(37) 

.(38) 


The scalar and gauge field equations of motion resulting from the Einstein-frame 
Lagrangian (here Lgff is considered function of (p,X,F‘^): 




\ (9Z/eff 


= 0 




dLes\ 


= 0 


(39) 


nv dx ) dip ^ dF^J 

thanks to the presence of the two (±) large flat regions (|33ll - (IM]l and ([36|) - 
as well as due to the “k-essence”-type nonlinear dependence of LeS on the scalar 
kinetic term, allow for the following two classes of nontrivial “vacuum” solutions: 

• (i) “Standard vacuum” containing standard constant “dilaton” vacuum plus 
nontrivial gauge field vacuum: 


p = const 


V = 0 


dLes 


= 0 


dLeS 


= 0 . 


(40) 


dp ’ 5F2 

Here the value p = const belongs to either the (—) flat region (IMll or the 
(+) flat region (l36ll of the effective scalar potential. 

• (ii) “Kinetic vacuum” (this type of “vacuum” exists thanks to the nonlinear 
w.r.t. X “k-essence” nature of the effective Lagrangian (l27ll l: 


dLes 

dX 


= 0 


dLeS 

dp 


= 0 


dLeS 

dF^ 


= 0 , 


(41) 
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In the first class of “standard vacuum” solutions the last equation (HOI) yields 
the following non-trivial “vacuum” value for the gauge field: 

V-Flc = , (42) 


and for the associated matter energy-momentum tensor (cf. (I22p ') we get: 

rpeS _ - r I _ rr(standard) 

^ HU — 9tJ,v^eS f!Letf =0“ 

’ aF 2 


(43) 


where total effective scalar potential in the “standard vacuum” 

dMl): 


Tj^(standard) _ rr i ^ 2 _ rr , ^ /o (^ 4eX2t4ff) 


In both (=f) flat regions (I001) - (I05]) and (1001) - 


and: 


J-Fl, = 


we have correspondingly: 

eVo 




^ ^ X 2 (1 + e/iV /2 + e^e/o) ’ 

eVo 


(+) X2 (1 + eMl/M2 + e^efi) ’ 




(standard) 


total 


~ u, 


(standard) _ 

(-) 


1 


1 - 


1 


fy'(standard) ^y(standard) 1 

^total — ^(+) 


4ex2 L" 1 + e/iV /2 + ee2/o 
1 


1 - 


4ex2L l + eM2/M2 + ee2/o"J ■ 


(44) 

(45) 

(46) 

(47) 

(48) 


Therefore, according to (|43ll the solutions of the Einstein-frame ( 7 ^, 4 -equations are 
of de Sitter type (pure de Sitter or Schwarzschild-de Sitter): 

\ 

ds^ = g^^dx^dx'^ = —A{r)dt^ + ^ + r'^i^dO^ + sin^ 9d(j)j , (49) 

A{r) = 1 - , or A{r) = 1- — - (50) 

3 r 3 

in static spherically symmetric coordinate chart, with effective cosmological con¬ 
stants A(±) given by (l47ll - (l48)) : 


A(_) = = 1 ^(^tandard) dSO]) - (l35| , (51) 

A(+) = = ic/(^t“dard) dMl) - jOHl) . (52) 


Erom the above analysis of the “standard vacuum” solutions - the one corre¬ 
sponding to ip = const belonging to the (—) fiat region of the effective scalar poten¬ 
tial with nonzero gauge field vacuum value and vacuum energy density as in (|45L 
dlZl), and the second the one corresponding to (p = const belonging to the (-I-) fiat 
region of the effective scalar potential with gauge field vacuum value and vacuum 
energy density as in (l46ll , (l4^ - we conclude that these “standard vacuum” solutions 
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describe charge confining phases with dynamically generated cosmological constants 

m, (ESI). 

Indeed, according to ‘t Hooft’s confinement proposal and as shown ex¬ 

plicitly in Ref. |28] in the case of flat spacetime, the non-zero vacuum values of 
the gauge field (l45l) . (l46l) imply confinement dynamics of charged particles with 
the strength of confinement proportional to these same gauge field vacuum values. 
Namely, under plausible truncation for static spherically symmetric configurations 
the canonically quantized theory in flat spacetime of charged fermions interacting 
with the nonlinear gauge fields with the “square-root” Maxwell term produces an 
effective “Cornell”-type potential I4ff = Eg. (1551) below) between 

quantized fermions separated by a distance L and where fi is proportional to the 
coupling constant of the “square-root” Maxwell term, which in turn is proportional 
to the non-zero vacuum value of the gauge field. 

The formalism to prove confinement used in Ref. [55] can be easily generalized 
to the case of curved static spherically symmetric spacetimes, in particular for de 
Sitter spacetime where both charged fermions are located within the interior de 
Sitter region below the de Sitter horizon (r < \/3/A(±)) - see Appendix B. 


4. “Kinetic Vacuum” Solutions 


We now turn our attention to the second class of “kinetic vacuum” solutions m- 
The equations = 0 and ^0- = 0 yield: 


Vkin — ~ 


“-^kin — ^eff 


1 - ex2/o/effe^ff 

2B 1 - ’ 

/eff - eX 2 foA‘^/B 




(53) 

(54) 


Using the identity = 2BX+A(l—ex2foV—F‘^) we can rewrite the Einstein- 
frame Lagrangian Lgg (12311 in the form: 


with: 


UefF — 


1 / dLes Y 

AB\ dX ) 




(55) 


U(v.,F2) = t/eff + ^-iv^(/eff-eX2/o^)-iF^(^-e"xl/o"^) • (56) 

Inserting in (I55]) - (I551) the “on-shell” values (|55)) - (IM1) . we obtain for the matter 
energy-momentum tensor: 

— 9livLeS |fLL!tf_g -0= ^ I (5"^) 

dX ’ dF'^ 

where total effective scalar potential in the “kinetic vacuum” (HTl) : 


Tr(kinetic) 

^total 


UeB 


4B 


^ 2 
reeff 


(/eff - eX2/o^) 

1 _ p2 ,2^2 f2A^ 
^ ®eff® X 2/0 B 


(58) 
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Note from Eos. (1551) - (I5S1) that within the “kinetic vacuum” = 0 the effective 
gauge coupling constants become: 

7 eff = /eff - eX2/o^ , efff = -- 2 2 f 2 


4.1. “Kinetic Vacuum” in the (+) Flat Region of Effective Scalar 
Potential 


Eirst we consider the “kinetic vacuum” solution in the (+) flat region. Using ((36ll - 
(l38)) in (15^ . (I53ll and (|58ll we obtain from (|59|) : 


/eff = /(+) = /(+) - 6X2/0 




B, 


which yields: 


_p 2 

1,; 


kin |(+) 


= 0 


Nlkin — -^(+) = ~ 


(+) 

A 


= 0 , 


(+) 


1 


2B. 


(+) 




(kinetic) 


total 


~ u, 


(kinetic) 


1 


jreff _ 


2 eX 2 

1 


(+) 4e^2 ' 

i.e., we have here an effective cosmological constant: 


4ex2 ’ 


A(+) = A 


(kinetic) 

(+) 


1 


8ex2 


(60) 


(61) 

(62) 


(63) 


Let us particularly stress on the first relation in (1611) - the zero vacuum value 
for the nonlinear gauge field, which is due to the vanishing (l60l) of the effective 
coupling constant of the “squ are-roo t” Maxwell term. Again, in accordance with ‘t 
Hooft’s confinement proposal 1 ^2 1 33 1 g^g demonstrated explicitly in Ref. [28] and 
in Appendix B the latter implies absence of confinement of charged particles, i.e., 
the “kinetic vacuum” (EU-dSSI) describes a deconfinement phase. 

According to (1571) and (l62ll - (l63ll the solutions of the Einstein-frame g;_ji/-equations 
in the “kinetic vacuum” are again of de Sitter type (l49ll - (150)) with A(_|_) given by 

(ESI). 

The equation for the “dilaton” “kinetic vacuum” (second Ea. (15T1) ') reads explic¬ 
itly: 


g^’'d^ipd^ip -= 0 . (64) 

6 X 2 

It has precisely the form of Hamilton-Jacobi equation for the Hamilton-Jacobi ac¬ 
tion: 


1 p^out I dr^ 

corresponding to spacelike geodesics x^{X) starting from some fixed point X(o) (e.g., 
X(o) = 0) at a fixed value of the affine parameter Ain and passing through x = a:(Aout) 
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at Aout- This Hamilton-Jacobi action (1651) measures the proper distance between the 
points a:(o) and x on the manifold modulo the numerical factor l/^ex 2 (for modern 
pedagogical exposition, see, e.g.. Ref. ED. 

A static spherically symmetric solution for ip{x) is given by: 




ex2Al(r 


(f{r) = (/?(+) 


dr' 


y/eX2 J y^A(r') 


( 66 ) 


where the initial value must belong to the (+) flat region (I551) - (I551) (large 
positive (f). 

In the case of pure de Sitter metric (l4^ - (l50ll the solution (IMp . measuring 
the proper radial distance between 0 and r, is clearly defined only for r in the 
interval r G (0,r(_|_)), where: 


i’(-H) = \/24ex2 


(67) 


is the de Sitter horizon radius. The solution <f{r) reads explicitly: 

ip(r) = ipi<\ + •\/^arcsin(^—) , (68) 

^(+) 

where the initial value belongs to the (j-) flat region of the effective scalar 
potential. Let us also recall that the integral in the second Ea. (l66l) . which is equal 
to arcsin(^:^), yields the proper radial distance in the internal de Sitter region 
r < r(+). 

In the case of Schwarzschild-de Sitter metric (H^ - dSOl) the solution i^{r) (1551) is 
defined in the interval r G {rs,rH) between the inner (Schwarzschild-type) horizon 
rs and the outer (de Sitter-type) horizon rgg. In what follows we will consider the 
case of pure de Sitter metric. 

Since the “kinetic vacuum” corresponding to the (-I-) flat region described by 
(|6TD-([68D is defined only within the finite-volume space region below the de Sitter 
horizon, in order to be extended to the whole space it must be matched to another 
spherically symmetric configuration with the standard constant “dilaton” vacuum 
defined in the outer region beyond the de Sitter horizon with: 

ip = v?(r(+)) = (p(+) J- V^TT = const for r > r(+) , (69) 

where the latter is the limiting value of (15^ at the horizon. The corresponding 
construction yields a gravitational bag-like solu tion mimicking both some of the 
features of the MIT bags in QCD phenomenologyas well as some of the features 
of the “constituent quark” model of Ref. |37] to be discussed in the next Section. 


4.2. “Kinetic Vacuum” in the (—) Flat Region of Effective Scalar 
Potential 

Next, let us consider the “kinetic vacuum” solution corresponding to the (—) flat 
region. In this case using (l33ll - (l35]l in (l54ll . (l53ll and (l5^ and introducing short¬ 
hand notations for some combinations of the parameters to simplify the resulting 
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[/_ 



Fig. 3. Qualitative shape of the “kinetic vacuum” energy density as function of ^ 

cs. 


expressions: 




7 = 


/2 

4e/f 


/3 = e 


f2 f2 

2 Jo J2 


2 -c2 2 


we obtain: 



■ e;^2(l + e-e^7)[l+^-e27(l + /3/7^)] ’ 

2 ex 2 l + C-e^7(l + W) ’ 

Ty'(kinetic) Ty-(kinetic) 1 ^ C C /^/T 

total - (-) = ^ l + e-e"7(l+OT) ’ 


(70) 

(71) 

(72) 

(73) 

(74) 


The spacetime metric is again of de Sitter type (l4^ - (l50)) with effective cosmological 
constant: 


A(_)^A 


(kinetic) 

(-) 


1 1 +^- eph 

8eX2 1 + ?-C^7(1 + /3/7^) ’ 


(75) 


According to Eq. (ITT]) the vacuum value of the gauge field is non-zero (except in 
the speci al case ^ = 0, see below (IMll b therefore, following ‘t Hooft’s confinement 
proposal I 32 | 33 | ^ Refs. [27]l3T] we conclude that the “kinetic vacuum” (ITTll - dT^ 

supports dynamics of charged particles as in the “standard vacuum” phase (l40ll 
except in the special case ^ = 0. 

The qualitative shape of the “kinetic vacuum” energy density in the (—) flat 
region (l7^ as function of the parameter ^ (1701) is depicted on 

Eig.3. 
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j^^(kinetic)has a local minimum at ^ = 0 , i.e., at 6 = 0 where: 

C/p“‘'c)(o) = i/4ex2 , (76) 

and it raises to +oo at ^ which are the roots of the quadratic expression 

1 +^-^^ 7(1 + / 3 / 7 ^: 



1 

27(1 + ^772) 


1 ± Vi +47(1+ ^772) 


(77) 


Alternatively, for large \^\, ± 00 ) = 4^(1 + l7ee^/o) ■ 

The scalar field for static spherically symmetric configurations is given by: 



(78) 


therefore, solutions exist only for the range of parameters for which A(_) < 0. From 
the explicit expression dZH) we find: 


X(_) < 0 for ^ ^ and for ^ < —2 , 


(79) 


where are the same as in (1771) . Thus, the solution of (1751) is similar to (1551) : 

/- r dr' ! - 

ip{r) = 93 (_) + y'2 |A(_)| J = ^i-) + v2|A(_)|r(_) arcsin(r7r(_)) , (80) 

where the initial value (/?(_) must belong to the (—) flat region (l33l) - (l35l) (large 
negative ip), is the de Sitter horizon radius: 


r(-) = Y^37A(_) (A(_)as in , (81) 

|A(_) I is given by (1751) . and again p{r) is defined only in the space region inside the 
de Sitter horizon (r < r(_)). 

Thus, similarly to the previous case for the (+) flat region, here in the (—) flat 
region again we need to match the “kinetic vacuum” given by dH])-® to another 
spherically symmetric configuration with the standard constant dilaton vacuum: 


p = ¥j(r(_)) 


P(-) + 7 rr(_) y^|A (_)|72 = const for r > r(_) 


(82) 


defined in the outer region, where p{ri^-)) is the limiting value of p{r) (I80|) at the 
horizon. This will be considered in the next Section. 

Let us specifically note that in the special case ^ = 0, i.e., for b = 

0 meaning that in this case the non-canonical scalar kinetic term is absent from 
the original second Lagrangian L(2) g])^ the expressions (TTn - dTSl) and (1751) . (15T1) 
drastically simplify: 


(7^1,-,) 


^(-) lj=o 2ex2 

* (kinetic) | _ 1 

Sex2 


Tr(kinetic) | _ ^ 

(-) l?=o 

l{=0= \/24eX2 .(83) 
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Expressions (IM|) precisely coincide with the corresponding values of |(+) = 

0, Jf(+), r(_|_) (l6T]) - (l63l) . (l67l) in the (+) flat region of the effective 

scalar potential. In particular, from the first relation in (1831) we conclude that in 
the special case ^ = 0 (6 = 0) the corresponding “kinetic vacuum” on the (—) flat 
region of the effective scalar potential implies deconfinement in complete analogy 
with the “kinetic vacuum” on the (+) flat region (subsection 4.1). 


5. Gravitational Bag-like Solutions 

5.1. Matching “Kinetic Vacuum” to Standard Constant “Dilaton” 
Vacuum in {+) flat region 

First, we construct matching of the “kinetic vacuum” in (+) flat region of the 
effective scalar potential given by de Sitter metric (l4^ - (l5n]) in the interior region 
(r < r(_|_)) below the de Sitter horizon r(_|_) = v^24ex2 with effective cosmological 
constant (|63|) and by Eqs. dSU)-®, to a static spherically symmetric configuration 
containing the standard constant “dilaton” vacuum (IMI) in the outer region (r > 
r(+)) beyond the de Sitter horizon. The “matching” specifically means that the 
“dilaton” field, the gauge field strength and the metric with its first derivatives 
must be continuous across the horizon, in particular, the de Sitter horizon of the 
interior metric must coincide with a horizon of the exterior metric. 

Obviously, the static spherically symmetric configuration in the outer region 
cannot be the “standard vacuum” of the full “dilaton” plus gauge field subsystem 
given by (|i(I)) . and (H51) . since: 

(a) In the inner “kinetic vacuum” region y — I (+)~ (first Ea. dCTl l. whereas 

in the outer “standard vacuum” region 7 ^ 0 (l46l) . which would imply that 

there should be a lightlike (“null”) brane with a non-zero surface electric charge 
density located on de Sitter horizon to account for the jump of the gauge field 
strength across the horizon; 

(b) The effective cosmological constant in the outer “standard vacuum” region 
(l5^ is different and smaller than the effective cosmological constant (l63ll in the 
inner “kinetic vacuum” region. 

In Refs. [42l |43] (see also the earlier works .^^0 have already explicitly de¬ 

rived static spherically symmetric solutions of the coupled gravity/nonlinear gauge 
field/scalar “dilaton” system (l27ll with a generalized Reissner-Nordstr 6 m-(anti)de 
Sitter geometry carrying a non-vanishing background constant radial electric field 
in addition to the standard Coulomb field. We will use this type of solution in the 
outer region beyond the de Sitter horizon to be matched with the “kinetic vacuum” 
(l 6 T|) - (| 68 |) in the interior region. r™™ 

Specifically, for r > r(+) = ^/24ex2 the solution reads (the additional 

factor in Ea. (l85|) below is due to the adopted normalization for the Newton 
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constant GNewton = l/167r): 


ds‘^ = -Aout{r)dt^ + 


dr^ 


■^out (l*) 


1 r 9m 

Aut(r) = 1 + —-V8^|Q|/(+)-h ^ ^ 

loTT L ^ ^ T ef 


+ sin^ 0d(j^ , (84) 
,(85) 


"(+)' 


t^out - I 


(standard) 


1 


1 - 


1 


lout - 2 ^ (+) 

\/-FLt{r) = V2\Eout{r)\ = e(+)/(+) - 


8ex2 L" 1 + eMf/M2 + ee2/2J 

\Q\ 


(as in (|48)) . (l5^ i 


\/^7 


(e(+) , /(+) as in ([381)) ,(86) 


= (/j(r(_|_)) = const (as in (IG^ I .(87) 


Notice that in (1551) we have taken the constant radial background electric field and 
the Coulomb field with opposite directions. Let us also point out that the scalar 
potential corresponding to the static radial electric field in (1551) : 

Elnt = -For = drAoir) , Ao{r) = -^e^+)/(+) r + , (88) 

which is a static spherically symmetric solution of the nonlinear gauge field equations 
(last Ea.(l39l) with as in (l27l) and with X — Q ~ “standard dilaton vacuum”), 
resembles the form of the well-known phenomenological “Cornell potential” in QCD 
which contains both a linear confining and a standard Coulomb part SlHini 

The mass and electric charge parameters (m, Q) in (l85ll are to be determined 
from the matching at the common horizon: 


>.^out (^(-t-)) — -^in)^)-!-)) — 0 , ^r^out (^(-t-)) — ^r'^in (^(-f-)) , (89) 

\l-Font{r{+)) = \J-F^i^ l(-H)= 0 (according to first Eq.dM])) , (90) 

with the “kinetic vacuum” dm-dMi) for r < r(_|_) = v^24ex2 with: 

dv^ \ 

ds^ = —Ain{r)dt^ + ——^ -\-r^{^d9^ + sin^ 9dcj)] , (91) 

A„(r) = l-^r2 ^ = ![/{“) = l/8ex2 . (92) 

Inserting (1551) into (1551) we determine Q: 

101 = 24.„ = i + ■ (93) 

Now, inserting (l93l) into (1551) yields: 

TO = 0 , (94) 


and the following relation between the integration constants Mi ^2 and the initial 
coupling constants e,e, fo- 


To recapitulate, we have obtained the following “vacuum-like” solution: 


(95) 
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• In the inner space region r < = y/24ex2 we have an interior de Sitter 

region eiD-dSSD below the de Sitter horizon at r = with effective cos¬ 
mological constant (|Mll . with vanishing vacuum gauge field (first Ea. (|6T]l i. 
“kinetic vacuum” scalar “dilaton” according to (|M|l and vacuum energy 
density (|62l) : 




1 

4ex2 ' 


(96) 


• In the outer space region r > r(_|_) = v'24ex2 we have static spherically 
symmetric metric (1851) with: 


-^out (j" ) — 1 


I 

2 ^ 


6 x 2 
^ e^Po 


1 


24ex2 


(^ 


I 

4^ 


) 


(97) 


where we have used the explicit expressions for /(+) (IMll . Q (IMll and 
^^(rtandard) together with the relation (IMl) . The metric with (|^ has 
de Sitter type horizon again at r = where the relation (|M]) among the 
parameters holds. 

• The outside nonlinear gauge field (l86l) is a static radial electric field of the 
explicit form: 


\/^M = ^K„.MI = j^(i-^). (98) 

where again we have used (IMll and (IMll . In (IMl) there is a Coulomb piece 
in addition to a non-zero background constant radial electric field: 


IE, 


background I 


= ^ef+)/(+) = 


■/24ex2/o 


(99) 


Thanks to the latter the Coulomb field is completely cancelled at the hori¬ 
zon. 

• The scalar “dilaton” is constant dlZl) and the energy density (p = —Tg) 
reads (using again (1551) . (H51) and (IMl) '): 


Pout(r-) 

46X2^^ 4ee2/2) 

Obviously (recall in (I100|) r > r(+) 



1 O 24ex2\ 
eeVo V 


= V24ex2): 


( 100 ) 


Poutir) < UP‘ 


(standard) 


4ex2 


1 ] 
l + eM2/M2 + ee2/o- 


< Pin = 


I 


4ex2 ' 
( 101 ) 


The above solution (IM1) - (I101I) is a gravitational bag-like configuration on the 
(-I-) flat region of th e effective scalar potential which mimics some of the properties 
of the MIT bagldSI Indeed, as already noticed in Section 3: 
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(i) In the inner finite volume space region below the horizon (r < the van¬ 
ishing vacuum value of t he g auge field (first Eg. (1611) 1 implies absence of confinement 
of charged particles I 221 - 131 J 

(ii) According to (IIOII) the vacuum energy density pin in the inner finite volume 

space region (for r < is larger than the energy density pout in the outside 

region. 

There are, however, other properties of the present gravitational “bag” solu¬ 
tion which are substantially different from those of the MIT bag and which rather 
resemble some of the properties of the solitonic “constituent quark” model 

(a) It is charged (the overall charge Q is non-zero (IMll l. 

(b) It carries non-zero “color” flux to infinite - because of the non-zero back¬ 
ground constant radial electric field (IMl) . 


5.2. Matching “Kinetic Vacuum” to Standard Constant “Dilaton” 
Vacuum in { — ) flat region 


Using the same procedure above we can construct the matching of the “kinetic vac¬ 
uum” (ffTll - dSnil in (—) flat region of the effective scalar potential given by de Sitter 
metric (l4^ - (|50)l in the corresponding interior region (r < r(_)) below the de Sitter 
horizon r(_) (1^ . with effective cosmological constant (fTfi]) . to a static spherically 
symmetric configuration containing the standard constant “dilaton” vacuum (I82() 
in the outer region (r > r(_)) beyond the de Sitter horizon. 

In the special case ^ = 0 (6 = 0), as already noted in (|83ll above, the expressions 
([7T1l-([73ll and (ESI) in the (—) fiat region precisely coincide with the corresponding 
expressions (|6T1l - (l63)l . (l67)l in the (-I-) fiat region. Taking into account (IMj) and the 
explicit form of e(_), /(_) (IM)l versus e(+), /(+) (IM]) and repeating the same steps as 
in subsection 5.1 we obtain in the special case 6 = 0 a completely analogous solution 
for the matching of interior “kinetic vacuum” dnli-dHQi) corresponding to the (— ) 
flat region of the effective scalar potential with the exterior region with the same 
generalized Reissner-Nordstrom-de Sitter geometry carrying a non-zero constant 
radial background electric field as in (IM1) - (I551) with (IM1) - (IM1) upon substitution 
(e(_i_),/(_!_)) —>■ (e(_),/(_)) and Mi _2 —t fi., 2 - In particular, instead of (1551) we now 
obtain the following relation among the parameters of the model: 

1 + - Se/^e^ = 0 . (102) 

h 

For the energy densities inequality instead of (IIOII) we now have: 


Pout(r) < 


4eX2 L 


1 - 


1 


< Pin = 


1 


4ex2 


(103) 


1 + e/?//2 + ee2/o J 

Thus, in the special case ^ = 0 (5 = 0) we obtain a gravitational bag-like configu¬ 
ration on the (—) fiat region of the effective scalar potential with the same properties 
as those of the gravitational bag solution in subsection 5.1 above - properties (i)-(ii) 
and (a)-(b). 
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One can straightforwardly extend the above solution to the general case of the 
parameter b 0 (i.e., ^ 7 ^ 0). However, the explicit expressions for the parameters 
Q and m in the exterior generalized Reissner-Nordstrom-de Sitter metric carrying 
a non-zero constant radial background electric field as well as the generalization 
of relation ( 11021 ) among the theory’s parameters become algebraically much more 
complicated. In particular, now the mass parameter m 7 ^ 0. Moreover, in the general 
case of 5 7 I 0 7 ^ 0 ) the vacuum value of the gauge field (1711) 1 in the inner finite 

volume space region below the de Sitter horizon (r < is non-vanishing, thus 

again implying confinement dynamics of charge particles as in the outer space region 
beyond the de Sitter horizon (r > r(_) (I 8 ip l where we have non-zero constant radial 
background electric field |^backgroundI = (the counterpart of ((M|) b On 

the other hand, according to (llO.jp the vacuum energy density pin in the inner finite 
volume space region is larger than the energy density pout in the outside region. 

Therefore, in the general case 6 7 ^ 0 7 ^ 0) the “vacuum-like” solution describing 

matching of “kinetic vacuum” to standard constant “dilaton” vacuum in the (—) 
flat region of the effective scalar potential is a gravitational bag-like solution which 
shares some of the properties of the “constituent quark” model (properties (a)-(c) 
in subsection 5.1 above), however, it does not at all mimic the properties (i)-(iii) of 
MIT bag unlike the gravitational bag-like solution of subsection 5.1. 

6. Conclusions 

In the present paper we have constructed a new kind of gravity-matter theory 
coupled to a nonstandard nonlinear gauge theory with the following noteworthy 
features: 

• Instead of the canonical Riemannian spacetime volume-form (generally co¬ 
variant integration measure density in terms of \/—g), to construct the 
action of the model we are employing two different and independent non- 
Riemannian spacetime volume-forms defined in terms of two auxiliary an¬ 
tisymmetric tensor gauge fields of maximal rank. 

• The action of our model contains apart from the standard Einstein-Hilbert 

R and Maxwell gauge field —F^ Lagrangian terms also scalar “dilaton” 
parts with a non-canonical kinetic term, as well as additional R? and a 
“square-root” Maxwell term The specific form of our action is dic¬ 

tated by the requirement of global Weyl-scale invariance. 

• Solving the equations of motion for the auxiliary antisymmetric tensor 
gauge fields building up the two non-Riemannian spacetime volume-forms 
introduces several arbitrary integration constants, some of them sponta¬ 
neously breaking the global Weyl-scale symmetry of the initial action. 

• The physical meaning of the above arbitrary integration constants is re¬ 
vealed within the canonical Hamiltonian formalism. Namely, these integra¬ 
tion constants turn out to be conserved Dirac-constrained canonical mo¬ 
menta conjugated to some of the components of the auxiliary antisymmetric 
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tensor gauge fields of maximal rank, the latter turning out to be essentially 
pure gauge non-propagating degrees of freedom. 

• After passing to the physical “Einstein frame” thanks to the appearance of 
the above arbitrary integration constants we obtain a remarkable effective 
matter Lagrangian of quadratic “k-essence” type. First, the latter contains 
an effective scalar “dilaton” potential of a very interesting form possessing 
two infinitely large flat regions for large negative and large positive “dila¬ 
ton” tp values. Second, all the remaining terms in the “k-essence” matter 
Lagrangian appear multiplied by nontrivial “dilaton”-dependent coefficient 
functions, including nontrivial effective gauge coupling constants running 
with If. 

• We study systematically the static spherically symmetric “vacuum”-like 
configurations corresponding to each of the flat regions of the effective scalar 
potential. 

• First, we find two globally existing in space phases corresponding to the 
standard constant “dilaton” vacuum values either in either of the two in¬ 
finitely large flat regions. These both phases describe confinement since in 
both cases the vacuum value of the nonlinear gauge field is non-zero. 

• Further, we find two localized in space (in “bubbles”) deconfining (confine¬ 
ment free) phases corresponding to the so called “kinetic dilaton vacuum” 
(when the quadratic “k-essence” effective action in the Einstein frame is 
extremized w.r.t. X - the “dilaton” kinetic term), where the vacuum value 
of the nonlinear gauge field vanishes. In one of these deconfining phases the 
“dilaton” lies on the (+) flat region of the effective scalar potential and in 
the second one it belongs to the (—) flat region with the additional restric¬ 
tion on the Lagrangian parameter b — 0. On the other hand, in the generic 
case of 6 ^ 0 the “kinetic dilaton vacuum” on the (—) flat region describes 
a localized confining phase since the vacuum value of the nonlinear gauge 
field is again non-zero there. 

• The localized deconfining phases inside the “bubbles” coexist with outside 
configurations corresponding to standard constant “dilaton” vacuums and 
nontrivial nonlinear gauge fields carrying non-zero constant radial back¬ 
ground electric field The energy density inside the “bubbles” is larger that 
the outside energy density. Thus, the full solution inside and outside the 
“bubbles” is a gravitation al ba g-like solution mimicking some of the basic 
properties of the MIT bag ^1^. 

• On the other hand, the analogy of the above gravitational bag with the 
MIT bag is only partial one. The present gravitational bag possess other 
properties (overall charge, carrying non-zero “color” flux to infinity) which 
resem ble some of the main properties of the solitonic “constituent quark” 
model 
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Appendix A. Canonical Hamiltonian Treatment of Gravity-Matter 
Theories with Non-Riemannian Volume-Forms 

Here we will briefly discuss the application of the canonical Hamiltonian formalism 
to the new gravity-matter model based on two non-Riemannian spacetime volume- 
forms O- In order to elucidate the proper physical meaning of the arbitrary inte¬ 
gration constants X 2 , Mi, M 2 o encountered within the Lagrangian formalism’s 
treatment of © it is sufficient to concentrate only on the canonical Hamiltonian 
structure related to the auxiliary maximal rank antisymmetric tensor gauge helds 
and their respective conjugate momenta. 

For convenience let us introduce the following short-hand notations for the field- 
strengths of the auxiliary 3-index antisymmetric gauge fields (the dot indi¬ 

cating time-derivative): 

$i(A) =A +a,A* , A = , A* = , (A.l) 

MB) =B+d,B^ , , B^ =, (A.2) 

<^{H)^H+dM , H = = , (A.3) 

Also we will use the short-hand notation: 

M\u,u) = RalM , L^^\u,u) = M + eR^ , (A.4) 

where are as in ©-© and where {u, u) collectively denote the set of the basic 
gravity-matter canonical variables (it) = A^j) and their respective velocities. 

For the pertinent canonical momenta conjugated to (IA.1I) - (IA.3I) we have: 

tta = Li{u,u) , ttb = L'''^\u,u) ++0^) , 

tth = -^{B +dM) , (A.5) 


TTAi =0 , TTgi =0 , TTHi = 0 . (A.6) 

The latter imply that A^,B^,H'^ will in fact appear as Lagrange multipliers for 
certain first-class Hamiltonian constraints (see Eas. (|A.10|) - (|A.lip below). Eor the 
canonical momenta conjugated to the basic gravity-matter canonical variables we 
have (using last relation (IA.5I) '): 

Pu = {A Zi(w,'u) + 7rH\/^—■ (A.7) 

o u an 

Now, relations (IA.5P and (IA.7I) allow us to obtain the velocities u, A, B, H 
as functions of the canonically conjugate momenta u=u (it,p„, ttb,tt//) etc. 
(modulo some Dirac constraints among the basic gravity-matter variables due to 
general coordinate and gauge invariances). Taking into account (IA.5p - (IA.6ll (and 
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the short-hand notations (IA.1|) - (|A.4I) I the canonical Hamiltonian corresponding to 

O: 


H =PuU +TrA A -t-TTs B +TrH H - {A +diA'')Li{u,u) 

+ —^{H +diW) 


(A.8) 


acquires the following form as function of the canonically conjugated variables (here 

'H=PuU 

+y/^TrHnB - diA'-TTA - diB^TVB - ■ (A.9) 


From (IA.9I) we deduce that indeed A'^,B'^, W are Lagrange multipliers for the first- 
class Hamiltonian constraints: 


diTTA = 0 —>■ tta = —Ml = const , (A. 10) 

and similarly: 

ttb = —M 2 = const , tth = X 2 = const , (A.11) 

which are the canonical Hamiltonian counterparts of Lagrangian constraint equa¬ 
tions of motion (jl4p . 

Thus, the canonical Hamiltonian treatment of o reveals the meaning of the aux¬ 
iliary 3-index antisymmetric tensor gauge fields A^^^x^ B^j^^x, H^^x - building blocks 
of the non-Riemannian spacetime volume-form formulation of the modified gravity- 
matter model CD - Namely, the canonical momenta tta, t^b, conjugated to the 
“magnetic” parts A, B, H (jA.ip - (IA.3p of the auxiliary 3-index antisymmetric tensor 
gauge fields are constrained through Dirac first-class constraints (lA.lOD - dA.llIl to 
be constants identified with the arbitrary integration constants X 2 , Mi, M 2 da 
arising within the Lagrangian formulation of the model. The canonical momenta 
■k\, conjugated to the “electric” parts A'‘,B'^,W (IA.ip - (|A.3p of the auxil¬ 

iary 3-index antisymmetric tensor gauge field are vanishing (IA.6I) which makes the 
latter canonical Lagrange multipliers for the above Dirac first-class constraints. 


Appendix B. “Cornell”-Type Confining Potential in Curved 
Spacetime 


Here we will follow the steps of the derivation in Ref. [28l of effective “Cornell”- 
type confining potential I44 | 45|46 | i^et^^egn quantized charged fermions based on the 
general formalism for quantization within the canonical Hamiltonian approach 
a’la Dirac of truncated gauge and gravity theories by imposing explicitly spherical 
symmetry on the pertinent Lagrangian action. 

In the present case the corresponding nonlinear gauge field action: 




L{F^) - 




fo 


. / F 2 
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yields equations of motion: 


, L'(f2) = --(i--^) , (B.2) 


whose fj, = 0 component - the nonlinear “Gauss law” constraint equation reads: 
1 






= 


V y/^J 


F' 


Oi 


(B.3) 


with D = (U®) denoting the electric displacement field nonlinearly related to the 
electric field E = (F°*) as in the last relation (IB.3p . 

In the special case of nonlinear gauge field theory (jB.ll) there exists a nontrivial 
vacuum solution y/—F^^^ = /g, which implies simultaneous vanishing of the electric 
displacement field, D = 0 meaning zero observed charge, and at the same time 
nontrivial electric field. In particular, for static spherically symmetric fields in static 
spherically symmetric metric (of the form (1491) with general A{r)) the only surviving 
component of is the nonvanishing radial component of the electric field £’’’ = 
—For, so that y/—f^ = v^lF] = /o- This can be viewed as the simplest classical 
manifestation of charge confinement: D = 0 and nontrivial E. 

Here we will employ the canonical Hamiltonian treatment in Ref. |47] and will 
truncate the nonlinear gauge field action to purely spherically symmetric fields, i.e., 
we will take For = d^Ar — drA^ independent of the space angles and the rest of the 
components of F^^, being zero. The action of the truncated theory reads: 


^truncated — 




drEnr^ 


\fI - A|Ff 


Or 


+ HoJ° +AJ" 


(B.4) 


L2*"” 

Note that in (IB.41) there is no explicit dependence on the Riemannian metric coeffi¬ 
cient A{r) = —goo = ilgrr- It is now straightforward to apply the canonical Hamil¬ 
tonian quantization procedure to (IB.41) within the Dirac formalism for constrained 
dynamical systems (e.g. Ref. H51). Obviously, in the case of de Sitter spacetime the 
radial coordinate r must be restricted to vary up to the de Sitter horizon radius rn- 
The canonically conjugated momenta w.r.t. Ao and Aj. read: 


H'J = 0 


(For - —^ 


V2J 


(B.5) 


where the first one H*^ = 0 is the standard primary Dirac constraint known in any 
gauge theory of Yang-Mills type. For the density of the canonical Hamiltonian one 
obtains: 

1 /o 


n = 


(H”)^ 


n'' + irr^fo - Arr + WdrAo - rA 


0 • 


(B.6) 


Sttc^ y /2 

Henceforth, for simplicity we will consider the case with no matter current J” = 0. 
Time-preservation of the primary constraint H^ = 0, i.e., ^H^ = = 0 

yields the standard secondary Dirac constraint - the “Gauss law” constraint: 

$i(r) = (9rn”-h J° = 0. (B.7) 
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Thus, one has to Dirac-canonically quantize the theory with canonical Hamilto¬ 
nian: 

1 /„^\2 


H = I dr 




(n-)' + An'+„yj 


(B.8) 


and with two first-class a’la Dirac constraints $o,i = 0 ($o = n° = 0 and = 0 
as in (jB.7() L which have to be supplemented by two canonically conjugate gauge¬ 
fixing conditions Xo,i- Since Aq and its conjugate momentum n° = 0 do not mix 
with the rest of the canonical variables they have no impact on the pertinent Dirac 
brackets between Ar and H'’ to be promoted to quantum operator commutators upon 
quantization. Thus we only need to choose an appropriate gauge fixing condition 
for the “Gauss law” constraint (jB.7(l . which we can take in the form: 


XI(r) = [ dz^Ax{z) . 
Jcir) 


(B.9) 


Here fc(r) path integral along a spacelike geodesic = x^(^) ending at the 
spacetime point with radial coordinate r. In particular, for the interior de Sitter 
region (r < rn) this spacelike geodesic can be taken in the form: 

TT 

= t = const , r(^) = rn sm^^/rn) , 0 < ^ < = r , 

(B.IO) 

where ^ is the de Sitter proper distance parameter, so that: 

pT 

Xi(r)=/ dzAr{z) , {«'i(r),xi(?’')| = 6(r - r') . (B.ll) 

Jo c P® 

Note that here and below 6{r — r') denotes the Dirac delta-function on the half-line 
(both r, r' > 0). 

It is now straightforward to calculate the Dirac bracket between the canonically 
conjugate pair given by: 

by using the standard Poisson bracket ^Ar{r), = 5{r — r'), which yields: 

|A,(r),n”(r')| =26ir-r'). (B.I3) 

Upon canonical quantization (IB.ljp becomes: 

n”(r),Hr(r')j = 2i<5(r-r') , i.e. U^r) =-2iS/SAr{r) . (B.I4) 

Now, following Ref. [28] we consider a gauge invariant quantum state of two op¬ 
positely charged (±eo) fermions located at r = 0 and r = L, respectively, explicitly 
given by: 

1$) = |'F(L)4'(0)) = ^(L)exp|ieo^ dzHr(z)|4'(0) |0) . (B.I5) 
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The average of the quantized canonical Hamiltonian (IB. 81) in this state (IB.15I) , where 
now n’’(r) will act on Ar{r) according to (IB.14I) : 

($|i7|$) = Feff(T) (B.16) 


can be viewed as effective potential between the quantized fermionic pair generated 
by the nonlinear gauge field theory containing the “square-root” Maxwell term 

dHU). 

Using (IB.141) one calculates: 


n’’(r),feo f dzAr{z) 

Jo 



dzAr(z) 


1 *eo 


dzAr(z) 


= 2eoO(L 
= SepiL 


r) , 
r) , 


(B.17) 

(B.18) 


where 9{r — r') denotes the step-function on the half-line (both r, r' >0). Plugging 
(IB.17I) - (IB.18I) into (IB.161) we obtain: 

g2 I 

V'eff(T) = + eo/oV^U + (L—independent const) , (B.19) 

Ztt E 


which has precisely the form of the “Cornell” potential 
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